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Abstract 



This paper presents an alternative pre-relativistic approach to the vacuum case of classical 
. electrodynamics represented by vacuum Maxwell equations. Our view is based on the under- 

standing that the corresponding differential equations should be dynamical in nature and the 
physical relations represented by them should be directly verifiable at least in principle, so they 
must represent local energy-momentum balance relations. 

One must consider the source of nonlinearity in a field theory 
' by starting with physically reasonable assumptions, rather 

than merely generalizing the purely mathematical form of the 
CN ! lagrangian or field equations. 

^ ; D.H.Delphenich 

o 

^ ; 1 Introduction: General Notions about Physical Objects and 

cn ! Interactions 
O 

In correspondence with modern theoretical view on classical fields we consider time dependent and 
space propagating electromagnetic fields as flows of spatially finite physical entities which have been 
called in the early 20th century photons. Since the efforts made during the past years to find ap- 
^ ■ propriate nonlinearizations of Maxwell vacuum equations [1], [2], [3], [4], [5], [6], [7], [8], [9], [10], [11], [12], 
, [13], [14], [15], and the seriously developed quantum theory have not resulted so far in appropriate, 

from our viewpoint, description of time stable entities of electromagnetic field nature with finite 
spatial support, in particular, adequate theoretical image of single photons as spatially finite physi- 
cal objects with internal dynamical structure, explaining in such a way their time stability and spin 
nature, we decided to look back to the rudiments of the electromagnetic theory trying to reconsider 
its assumptions in order to come to equations giving appropriate solutions, in particular, solutions 
with spatially finite carrier at every moment of their existence and space-propagating as a whole, 
keeping, of course, their physical identity. This reconsideration brought us to the necessity to ad- 
mit that introducing as clear as possible notions about the concepts physical object and physical 
interaction is strongly needed. Here we briefiy present the view on these concepts that we shall 
follow throughout the paper. 

When we speak about physical objects, e.g. classical particles, solid bodies, elementary particles, 
fields, etc., we always suppose that some definite properties of the object under consideration do not 



*e-mail:sdonev@inrne. bas.bg 



1 



change during its time-evolution under the influence of the existing environment. The availability 
of such time-stable features of any physical object makes it recognizable among the other physical 
objects, on one hand, and guarantees its proper identification during its existence in time, on the 
other hand. Without such an availability of constant in time properties (features), which are due to 
the object's resistance and surviving abilities, we could hardly speak about objects and knowledge at 
all. So, for example, two classical mass particles together with the associated to them gravitational 
fields survive under the mutual influence of their gravitational fields through changing their states 
of motion: change of state compensates the consequences of the violated dynamical equilibrium 
that each of the two particles had been established with the physical environment before the two 
gravitational fields have begun perturbating each other. 

The above view implies that three kinds of quantities will be necessary to describe as fully as 
possible the existence and the evolution of a given physical object: 

1. Proper (identifying) characteristics, i.e. quantities which do NOT change during the 
entire existence of the object. The availability of such quantities allows to distinguish a physical 
object among the other ones. 

2. Kinematical characteristics, i.e. quantities, which describe the allowed space-time evo- 
lution, where " allowed" means consistent with the constancy of the identifying characteristics. 

3. Dynamical characteristics, i.e. quantities which are functions (explicit or implicit) of the 
proper and of the kinematical characteristics. 

Some of the dynamical characteristics must have the following two important properties: they 
are in a definite degree universal, i.e. a class of physical objects (may be all physical objects) 
carry nonzero value of them (e.g. energy- momentum) , and they are conservative, i.e. they may 
just be transferred from one physical object to another (in various forms) but no loss is allowed. 

Hence, the evolution of a physical object subject to bearable/acceptable exterior influence 
(perturbation), coming from the existing environment, has three aspects: 

1. constancy of the proper (identifying) characteristics, 

2. allowed kinematical evolution, 

3. appropriate exchange of dynamical quantities with the physical environment guaranteing 
object's time stability. 

Moreover, if the physical object under study is space-extended (continuous) and demonstrates 
internal structure and dynamics, i.e. available interaction of time-stable subsystems, it should 
be described by a many-component mathematical object, e.g. vector valued differential form, 
clearly, we must consider this internal exchange of some dynamical characteristics among the various 
subsystems of the object as essential feature determining in a definite extent object's integral 
appearence. 

The above features suggest that the dynamical equations, describing locally the evolution of 
the object, may come from giving an explicit form of the quantities controlling the local internal 
and external exchange processes, i.e. writing down corresponding local balance equations. Hence, 
denoting the local quantities that describe the external exchange processes hy Qi,i = 1,2, ... , the 
object should be considered to be Qj-free, i = 1, 2, . . . , if the corresponding integral values are time 
constant, which can be achieved only if Qi obey differential equations presenting appropriately (im- 
plicitly or explicitly) corresponding local versions of the conservation laws (continuity equations). 
In case of absence of external exchange similar equations should describe corresponding internal 
exchange processes. The corresponding evolution in this latter case may be called proper evolution. 

Summerizing, we may assume the rule that the available changes of a physical field system, or a 
recognizable subsystem of a given system, must be refered somehow to the very system, or to some 
of the susbsystems, in order to evaluate their significance: 
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-if the refered quantity is zero, then the changes are admissible and the system/subsystem keeps 
its identity; 

-if mutually refered quantities among subsystems establish dynamical equilibrium, i.e., each 
subsystem gains as much as it loses, then the whole system keeps its identity; 

-if the refered quantity is not zero then the identity of the system is partially, or fully, lost, so, 
our system undergoes essential changes leading to becoming subsystem of another system, or to 
destruction, giving birth to new system(s). 

In trying to formalize these views it seems appropriate to give some initial explicit formulations 
of some most basic features (properties) of what we call physical object, which features would lead us 
to a, more or less, adequate theoretical notion of our intuitive notion of a physical object. Anyway, 
the following properties of the theoretical concept "physical object" we consider as necessary: 

1. It can be created during finite period(s) of time. 

2. It can be destroyed during finite period(s) of time. 

3. It occupies finite 3-volume at any moment of its existence, so it has spatial structure and 
may be considered as a system consisting of two or more interconnected subsystems. 

4. It has a definite stability to withstand definite external disturbances. 

5. It has definite conservation properties. 

6. It necessarily carries sufficiently universal measurable quantities, e.g. energy-momentum. 

7. It exists in an appropriate environment (called usually vacuum) , which provides all necessary 
existence needs. Figuratively speaking, every stable physical object lives in a dynamical equilibrium 
with the outside world, which dynamical equilibrium may be realized in various regimes. 

8. It can be detected by other physical objects through allowed exchanges of appropriate 
physical quantities, e.g. energy-momentum. 

9. It may combine/coexist through interaction with other appropriate physical objects to form 
new objects/systems of higher level structure. In doing this it may keep its identity and can be 
recognized and identified throughout the existence of the system as its constituent/subsystem. 

10. Its destruction gives necessarily birth to new objects, and this process respects definite rules 
of conservation. In particular, the available interaction energy among its subsystems may transform 
entirely, or partially, to kinetic one, and carried away by the newly created objects/systems. 

The property to be spatially finite we consider as a very essential one. So, the above features 
do NOT allow the classical material points and the infinite classical fields (e.g. plane waves) to 
be considered as physical objects since the former have no structure and cannot be destroyed at 
all, and the latter carry infinite energy, so they cannot be finite-time-created. Hence, the Born- 
Infeld "principle of finiteness" [2] stating that "a satisfactory theory should avoid letting physical 
quantities become infinite" may be strengthened as follows: 

All real physical objects are spatially finite entities and NO infinite values of the 
physical quantities carried by them are allowed. 

Clearly, together with the purely qualitative features, physical objects carry important quantita- 
tively described physical properties, and any external interaction may be considered as an exchange 
of the corresponding quantities provided both the object and the corresponding environment carry 
them. Hence, the more universal is a physical quantity the more useful for us it is, and this moment 
determines the exclusively important role of energy-momentum, which modern physics considers 
as the most universal one, i.e. it is more or less assumed that: 

All physical objects necessarily carry energy- momentum and most of them are able 
in a definite extent to lose and gain energy-momentum. 
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The above notes more or less say that we make use of the term "physical object" when we 
consider it from integral point of view, i.e. when its stability against external perturbations is 
guaranteed. We make use of the term "physical system" when time-stable interacting subsystems 
are possible to be recognized/identified and the behavior of the system as a whole, i.e. considered 
from outside, we try to consider as seriously dependent on its internal dynamical structure, i.e. on 
an available stable interaction of its time-recognizable/time-identifiable subsystems. Therefore we 
shall follow the rule: 

Physical recognizability of time-stable subsystems of a physical system requires 
corresponding mathematical recognizability in the theory. 

We note that further wc approach the problem from nonrelativistic viewpoint. 

Since any physical interaction presupposes equilibrium or nonequilibrium dynamical exchange, 
i.e. flows of some physical quantities among the subsystems, the above notes make us turn to the 
intrinsic dynamical nature of vector fields (on M^) as formal local representatives of real flows, or 
real dynamical equilibrium. 

2 Maxwell stress tensors generated by one and two vector fields 

Prom formal point of view the most natural mathematical object to be considered as stress gen- 
erating, seems to be any vector field, because every vector field defined on an arbitrary manifold 
M generates 1-parameter family (pt of (local in general) diffeomorphisms of M. Therefore, having 
defined a vector field X on M we can consider for each t G R the corresponding diffeomorphic image 
iftiU) of any region U C M. Hence, if the spatially flnite object under consideration occupies for 
each t the finite region Ut C M^, carries dynamical structure and propagates in the 3-space, should 
stay identical to itself during propagation, the available 1-parameter group of diffeomorphisms can 
be considered as the appropriate mathematical tool to describe its propagation as a whole, as well 
as, from local point of view, since every admissible change described by ipt has unique opposite, 
described by ip-t, i.e. nothing essential has been lost irreversibly and the object stays identical to 
itself. 

Let now X be a vector field on the euclidean space (M^,y), where g is the euclidean metric in 
TM^ having in the canonical global coordinates (x^, x^, = x, y, z) components gu = 522 = 533 = 1 
and 5^12 = gi3 = g23 = 0. The induced euclidean metric in T*]R^ has in the dual bases the same 
components and will be denoted further by the same letter g. The corresponding isomorphisms 
between the tangent and cotangent spaces and their tensor, exterior and symmetric products will 
be denoted by the same simbol g, so (summation on the repeating indecies is assumed) 

Having another vector field Y on we can form their scalar product 

g{X, Y) = X.Y = gijX'Y^ = X,Y' = X.Y' + XgF^ + x^Y^ 

According to classical vector analysis on [16] for the differential of the function g{X, Y) we can 
write 

dg{X, Y) = {X.V)Y + {Y.V)X + X XTot{Y) + Y x rot {X), 
where in our coordinates 

d dY^ d 

X.V = Vx=X'^, {X.V)Y = VxY = X'^^, 

ox^ ax* ox^ 
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" x" is the usual vector product, and 

{lotxy 

The Hodge *g-operator acts in these coordinates as follows: 



(dX^_dX^ dX^_dX^ dX^_dX^ 
\ dx^ dx^ ' dx^ dx^ ^ dx^ dx^ 



*dx = dy A dz, *dy = —dx A dz, *dz = dx A dy, 

*{dx A dy) = dz, *{dx A dz) = —dy, *{dy A dz) = dx, 
*{dx A dy A dz) = 1, *1 = dx Ady A dz. 
Corollary. The following relation holds (d denotes the exterior derivative): 

votX = (g)-^ * dg{X), or ^(rotX) = *dg{X). 

Assume now that in the above expression for dg{X, Y) we put X = Y. We obtain 

^dg{X,X) = ^d{X^) = X X rotX + {X.V)X = X x rotX + Vx^- 
In components, the last term on the right reads 

{VxXy = X'ViX^ = ViiX'X^) - xNiX' = Vi{X'X^) - XMlvX, 

where (denoting by Lx the Lie derivative along X) 

( dX^ \ dX^ 

divX = *Lx {dx Ady A dz) = * I -—^dx Ady Adz 



^ dx'^ J dx^ 

Substituting into the preceding relation, replacing d(X^) by (Vi6jX'^)dx^ and making some ele- 
mentary transformations we obtain 

Vi (^X'X^ - ^g'^X^^ = [(rotX) X X + XdivXy. 

The symmetric 2-tensor 

2 ^ 

we shall call further Maxwell stress tensor generated by the (arbitrary) vector field X G X(M'^). 
Clearly, when we raise and lower indices in canonical coordinates with g we shall have the following 
component relations: 

Mij = Mf = M'^ 

which does not mean, of course, that we equalize quantities being elements of different linear spaces. 

We note now the easily verified relation between the vector product " x" and the wedge product 
in the space of 1-forms on M^: 

XxY = ig)-^ (*(5(X)Aff(y))) 

= (g)-'^ o i{X A Y){dx Ady A dz), X,Y ^ X(M'^). 
We shall prove now the following 
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Proposition. If a = g{X) then the following relation holds {i{X)da means X^daijdx^): 

g{rot X X X) = i{X)da = - * {a A *da). 

Proof. 

g{mt X X X) = go * (^(rotX) A g{X)) = - * (a A *da). 

For the component of i{X)da before dx we obtain 

\dx^ dx^ J \dx^ dx^ J ' 

and the same quantity is easily obtained for the component of [— * {a A *da)] before dx. The same 
is true for the components of the two 1-forms before dy and dz. The proposition is proved. 

Hence, da = dg{X) is the 2-form across which the vector field X will drag the points of the 
finite region U C M^, and local interaction between X and da will take place if i{X)da is not zero. 

As for the second term XdivX of the divergence VjM^dx^, since d * a = divX{dx Ady A dz), 
we easily obtain 

i{g~^ {*a)d * a = X divX. 

Hence, additionally, the 2- vector g~^{*a) will drag the points of C/ C across the 3-form d * a. 
We can write now 

ViM'jdx^ = [i{X)da + i{g-^{*a))d * a].dx^. 

One formal suggestion that comes from this relation is that the interior product of a multivector 
and a differential form is appropriate quantity to be used as a quantitative measure of physical 
interaction. 

Another suggestion is that this relation may be understood in the sense, that our physical field, 
initially represented by the vector field X, can be equally well represented by two objects: a = g{X) 
and *a. Moreover, 

-there is no available local nonzero interaction stress between a and *a, 

-each of them interacts with itself only, so that there is no stress exchange between a and *a. 

The naturally isolated two terms in ViM^dx^ suggest: any realizable dynamical process repre- 
sented by the vector field X to be described by the two objects (a, *a), and the time-recognizable 
nature of a and *a to be guaranteed by the following equations: 

i{X)da = Q, z(5"^(*a))d*a = 0, i.e. X x rotX = 0, divX = 0. 

In the above consideration the role of the euclidean metric g was somehow overlooked since 
no derivatives of g appeared explicitly. We are going now to come to these equations paying due 
respect to g . 

We shall work further with the 1-form a = g{X), where the vector field X generates stress 
according to the Maxwell stress tensor Mj(X). We want to see how the 3-form a A *a = 
i{X)a.dx A dy A dz = i{X)a.uj changes along an arbitrary vector field Y ^ X, so we have to 
find the corresponding Lie derivative. 

Lyia A *a) = {Lyoi) A *a -|- a A Ly(*a) 

= {Lyoi) A *a + a A ^Lya + a A [Ly , *i]a = 2{Lyoi) A *a + a A [Ly, *i]a. 
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where [Ly , *i] is the commutetor Ly o *i — *i o Ly, and the index of * denotes the degree of the 
form it is apphed to. Further we get 



(Lya) A *q; = (iyda) A*a — {a, Y).d * a + d.{{a, Y). * a). 

Noting that Ly(a A *a) = d(Q:^iya;) and the relation between the exterior derivative d and the 
coderivative S in euclidean case for p-forms, given by 

(-l)P('^-p) Sp o = *(„_p+i) o d(„_p) 

we obtain consecutively 

d ^aid^iyu — {a, Y) * = — [a*(da)jj + aj divXjY^oj + A [Ly, *i]a, 

S 0*2 ( ^"jct* — (") Y) * a \ = —[a^{da)ij + aj divXjy-' + A [Ly, *i]a, 



*2(q!, y) * a = (a, Y).a, 



6 ( -aio'gjk - ajak)Y^dx 



aia^Sl - ako' ) Y 



So we get 



-aiO^Si - akoc' ] Y 



[a\da)ij + aj divXjy'' + *-a A [Ly, *i]a. 



Hence, if [Ly, *i] = 0, for example when K is a Killing vector field for g, and the equations for our 
1-form a = g{X) are 

a\da)ij = 0, aj divX = 0, i,j = 1,2,3 

then the co-closed 1-form MijY^dx^ defines the closed 2-form *{MijY'^dx^), so, its integral over a 
closed 2-surface that separates the 3-volume where, by supposition, a is diff'erent from zero, gives 
a conservative quantity and the nature of this conservative quantity is connected with the nature 
of the vector field Y, which is a Killing vector in our 3-dimensional case. Of course, the term 
"conservative" here is, more or less, trivial, since just stress is available, no space propagation takes 
place. 

Finally we note that if the vector fields Xi,X2, ■■■,Xp describe time-stable physical exchanges, 
we can form the sum of their Maxwell stress tensors 



P r 



M'^{Xi) + M'\X2) + ... + M'\Xp) = 



{X^nXkY - ^{Xkfg'^ 



k=l 

which is not equal to the Maxwell stress tensor of their sum M^^{Xi -|- ... -|- Xp). 
As for solutions of the above equations, we may suppose two classes of solutions: 

1. Linear, i.e. those satisfying da = 0, d*Q; = 0. These are generated by a function / satisfying 
the Laplace equation A/ = with a = df . 

2. Those, satisfying da ^ but a^daij = and d * a = 0. Note that for such solutions the 
determinant c/et||(da)ij||,z,_7 = 1,2,3, is always zero, so the homogenious linear system 



a\x,y,z)daij{x,y,z) = 
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allows to express explicitly X through the derivatives of its components. 

Of course, all these solutions are static, so they can not serve as models of spatially propagating 
finite physical objects with dynamical structure. 

We pass now to the case of two vector fields. 

Let V and W be two vector fields on our euclidean 3-space. Summing up the corresponding 
two Maxwell stress tensors we obtain the identity: 



= [(rotF) X y + l/divy + (rotW^) xW + WdwWy. 
Let now {a{x, y, z),b{x, y, z)) be two arbitrary functions on M^. We consider the transformation 

{V,W) {Va-Wb,Vb + Wa). 

Corollary. 

The tensor M(y^^^Y) transforms to (a^ + b^)M(y y/y 
Corollary. 

The transformations (F, W) — ^ {V a — W b,V b + W a) do not change the eigen directions 
structure of M^^y-^^y 
Corollary. 

If o = cos a,b = sin a, where a = a{x, y, z) then the tensor M(y ^y) stays invariant: 

^{V,W) = M {V COS a-W sin ay sin a+W COS a)- 

The expression inside the parenteses above, denoted by M*-^, looks formally the same as the 
introduced by Maxwell tensor M*-' (E, B) from physical considerations concerned with the elec- 
tromagnetic stress energy properties of continuous media in presence of external electromagnetic 
field (E, B). Hence, any vector or any couple of vectors (V,!^), defines such tensor which we 
denote by My, or M(y^yi/y and call Maxwell stress tensor. The term, "stress" in this general 
mathematical setting could be justified in the following way. Every vector field on generates 
corresponding flow by means of the trajectories started from some domain Ut=o C M^, where t is an 
arbitrary parameter: at t > the domain Ut=o is diffeomorphically transformed to a new domain 
Ut C M^. Having two appropriate vector fields on we may obtain two compatible flows, i.e. to 
expect the points of an interesting for us domain Ut=Q C to be forced to accordingly move to 
new positions inside a new appropriate domain. 

We emphasize the following moments: 

1. The identity satisfied by M(y^m^ is purely mathematical, and t in Ut is arbitrary parameter, 
not time in general; 

2. On the two sides of this identity stay well defined coordinate free quantities; 

3. The tensors M(yT^) do not introduce interaction stress: the full stress is the sum of the 
stresses generated by each one of the couple (F, W). 

Physically, we may say that the corresponding physical medium that occupies the spatial re- 
gion Uo and is parametrized by the points of the mathematical subrcgion Uq C M^, is subject to 
compatible and admissible physical "stresses" generated by physical interactions mathematically 
described by the vector fields iy,W), and these physical stresses are quantitatively described by 
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the corresponding physical interpretation of the tensor M^y^^^-^ . Clearly, we could extend the couple 
{V, W) to more vectors (Vi, V2, Vp), but then the mentioned invariance properties of M(y^^r•) may 
be lost, or appropriately extended. 

We note that the stress tensor M*-^ appears as been subject to the divergence operator, and if 
we interpret the components of M^^ as physical stresses, then its divergence acquires, in general, 
the physical interpretation of force density. Of course, in the static situation as it is given by the 
relation considered, no stress propagation is possible, so at every point the local forces mutually 
compensate: V^M*^ = 0. In electrodynamics, if propagation is allowed, then the force field may 
NOT be zero: VjM*-' ^ 0, and we may identify it as a real time-change of appropriately defined 
momentum density P(e,b)- So, assuming some expression for this momentum density P(e,b) are 
ready to write down corresponding dynamica field equation through equalizing the spatially directed 
force densities VjM*^ with the momentum density changes along the time variable, i.e., equalizing 
ViM*-' with the (rf)-derivative of P(e,b)) where c = const is the translational propagation velocity 
of the momentum density flow of the physical system considered. In order to flnd how to choose 
P(E,B) ™ case of free EM-field we have to turn to the intrinsic physical properties of the field. As 
the past years showed, namely is the appropriate carrier of the physical properties of the 

field, so, it seems natural to turn to the eigen properties of -A<^(eb)" 

3 Eigen properties of Maxwell stress tensor 

We consider M*-'(E,B) at some point p G and assume that in general the vector fields E and 
B are lineary independent, so E x B 7^ 0. Let the coordinate system be chosen such that the 
coordinate plane {x,y) to coincide with the plane defined by E(p),B(p). In this coordinate system 
E = {El, E2,0) and B = {Bi,B2, 0), so, identifying the contravariant and covariant indices through 
the Euclidean metric 5^^ (so that M*-^ = Mj = Mjj), we obtain the following nonzero components 
of the stress tensor: 

Ml = {E^f + {B^f - ^(E^ + B^); = Ml = E^ E2 + Bi B'^; 

M| = {E'^f + {B'^f - ^(E^ + B^); M| = -^(E^ + B^). 

Since = -M|, the trace of M is Tr(M) = -^(E^ + B^). 
The eigen value equation acquires the simple form 

[iMlf-i\f]+iM^f]{Ml-X) = 0. 

The corresponding eigen values are 

Ai = -1(E2+B2); A2,3 = ±^J{Mly + (M^y = ±^ VUi)^ + (^2)^ 

where Ii = B^ - E^, I2 = 2E.B. 

The corresponding to Ai eigen vector Zi must satisfy the equation 

E(E.Zi) + B(B.Zi) = 0, 

and since (E,B) are lineary independent, the two coefficients (E.Zi) and (B.Zi) must be equal to 
zero, therefore, .^1 7^ must be orthogonal to E and B, i.e. Zi must be colinear to E x B: 



9 



The other two eigen vectors Z^fi satisfy correspondingly the equations 



E(E.Z2,3) + B(B.Z2,3) = [ ± \y/{hf + {hf + ^(E' + B2)] ^2,3. (*) 



Taking into account the easily verified relation 



( 



E^ + B- 
2 




- IE X Bp 



so that 



E2 + B2 
2 



- |E X B| > , 



wc conclude that the coefficient before ^2,3 on the right is always different from zero, therefore, 
the eigen vectors Z2,3(p) lie in the plane defined by (E(p),B(p)), jy G M'^. In particular, the above 
mentioned transformation properties of the Maxwell stress tensor M{V, W) — ?• (a^ + b'^)M(V, W) 
show that the corresponding eigen directions do not change under the transformation {V, W) — >■ 
{Va-Wb,Vb + Wa). 

The above consideration suggests: the intrinsically allowed dynamical abilities of the field are: 
translational along (E x B) . and rotational inside the plane defined by (E, B), hence, we may expect 
finding field objects the propagation of which shows intrinsic local consistency between rotation and 
translation. 

It is natural to ask now under what conditions the very E and B may be eigen vectors of 
M(E, B)? Assuming A2 = + {I2Y and Z2 = E in the above relation (*) and having in 

view that E x B 7^ we obtain that ^(F?) + B(E.B) must be proportional to E, so, E.B = 0, i.e. 
I2 = 0. Moreover, substituting now /2 = in that same relation we obtain 



The case "-" sign before the square root, i.e. A3 = —\\/ {IiY + (-^2)^, leads to analogical conclusions 
just the role of E and B is exchanged. 

Corollary. E and B may be eigen vectors of M(E, B) only if Ii = I2 = 0. 

The above notices suggest to consider in a more detail the case A2 = — A3 = for the vacuum 
case. We shall show, making use of the Lorentz transformation in 3-dimensional form that, if 
these two relations do not hold then under E x B 7^ the translational velocity of propagation is 
less then the speed of light in vacuum c. Recall first the transformation laws of the electric and 
magnetic vectors under Lorentz transformation defined by the 3- velocity vector v and corresponding 
parameter /3 = v/c,v = \-v\. If 7 denotes the factor l/y^l — /3'^ then we have 



Assume first that I2 = 2E.B = 0, i.e. E and B are orthogonal, so, in general, in some coordinate 
system we shall have E x B 7^ . 

If Ii > 0, i.e. |E| < |B|, we shall show that the conditions E' = 0,v.B = 0, 00 > 7 > 
are compatible. In fact, these assumptions lead to 7V.E + (1 — 7)(E.v) = 0, i.e. E.v = 0. 




E' = 7 E H 2-^v(E.v) + -V X B, 



B' = 7 B + — 2"^v(B.v) - -V X E. 
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Thus, c|E| = v|B||sin(v, B)|, and since v.B = then |sin(v, B)| = 1. It follows that the speed 

lEI 

V = cjgj < c is allowed. 

If Ii < 0, i.e. |E| > |B|, then the conditions B' = and v.E = analogically lead to the 

IB! 

conclusion that the speed v = cj^ < c is allowed. 

Assume now that I2 = 2E.B 7^ 0. We are looking for a reference frame K' such that E' x B' = 0, 
while in the reference frame K we have E x B 7^ 0. We choose the relative velocity v such that 
v.E = v.B = 0. Under these conditions the equation E' x B' = reduces to 

E X B + -(E^ + B2) = 0, so, - = IE X B|/(E2 + b2). 
c c 

Now, from the above mentioned inequality E^ + B^ — 2|E X B| > it follows that f < 1- 

Physically, these considerations show that under nonzero Ii and I2 the translational velocity 
of propagation of the field, and of the stress field energy density of course, can be eliminated by 
passing to appropriate frame of mass bodies. Hence, the only realistic choice for the vacuum case, 
where this velocity is assumed by definition to be equal to the vacuum light speed c and can 
not be eliminated by passing to appropriate physical frame, is Ii = I2 = 0, which is equivalent 
to E^ + B^ = 2|E X B|. In view of this and assuming \Tr{M)\ = 5(E2 + b2) to be the stress 
energy density of the field, the names "electromagnetic energy flux" for the quantity cE x B, and 
"momentum" for the quantity ^E x B, seem well justified without turning to any dynamical field 
equations. 

These considerations suggest also that if /i = 0, i.e. |Ep = |Bp during propagation, then the 
energy density can be presented in terms of each of the two constituents, moreover, in this respect, 
both constituents have the same rights. Therefore, a local mutual energy exchange between the 
two subsystems formally represented in terms of E.B is not forbidden in general, but, if it takes 
place, it must be simultaneous and in equal quantities. Hence, under zero invariants /i = and 
I2 = 2E.B = 0, internal energy redistribution between the two subsystems of the field would be 
allowed but such an exchange should occur without available interaction energy because the full 
energy density is always equal to the sum of the energy densities carried by the two subsystems of the 
field. But the above relation E^ + B^ = 2|E x B| and the required time-stability and propagation 
with velocity "c" of the field suggest /imply also available internal momentum exchange to take 
place, since this relation means that the energy density is always strongly proportional to the 
momentum density magnitude ^|E x B|. However, the important observation here is that, verious 
combinations constructed out of the constituents E and B, e.g., {'Ecosa~'Bsina,Esina + 'Bcosa), 
where a{x, y, z, t) is a functon, may be considered as possible representatives of the two recognizable 
subsystems since they carry the same energy |(E^ + B^) and momentum ^|E x B|. Therefore, 
the assumption that the very E and B may be considered as appropriate mathematical images 
of recognizable time-stable subsystems of a time-dependent electromagnetic field does NOT seem 
adequate and has to be reconsidered. 

Hence, which combinations of E and B deserve represent mathematically the two subsystems 
of a time-dependent and space-propagating electromagnetic (EM) field? 

In view of these considerations we assume the following understanding: 

Every real EM-field is built of two recognizable subsystems, the mathematical 
images of which are not the very (E, B), but are expressed in terms of (E, B), both 
these subsystems carry always the same quantity of energy-momentum, guaranteeing 
in this way that the supposed internal energy- momentum exchange will also be in 
equal quantities and simultanious. 
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4 Real EM-fields viewed as built of two recognizable and perma- 
nently interacting subsystems. 

In accordance with the above assumption the description of dynamical and space-propagating 
behavior of the field will need two appropriate mathematical objects to be constructed out of the 
two constituents (E, B). These two mathematical objects must meet the required property that 
the two physical subsystems of the field carry always the same quantity of energy-momentum, 
and that any possible internal energy-momentum exchange between the two subsystems shall be 
simultaneous and in equal quantities. 

We are going to consider time dependent fields, and begin with noting once again that the 
assumption that the energy density of the field coincides with |Tr(M)| = ^[E^ + B^] suggests 
that there is NO interaction energy between the two subsystems of the field: the full stress tensor 
(and the energy density, in particular) is a sum of the stress tensors carried separately the two 
subsystems. As we mentiond above, this does NOT mean that there is no energy exchange between 
the two subsystems of the field. 

Now, following the above stated idea we have to find two appropriate mathematical images of 
the field which images are NOT represented directly by the electric E and magnetic B vectors, 
but are constructed out of them. In terms of these two appropriate mathematical representatives 
of the corresponding two partnering subsystems we must express the mentioned special kind of 
energy-momentum exchange, respecting in this way the fact that neither of the two constituents is 
able to carry momentum separately. 

In view of the above we have to assume that the field keeps its identity through adopting 
some special and appropriate dynamical behavior according to its intrinsic capabilities. Since the 
corresponding dynamical/field behavior must be consistent with the properties of the intrinsic 
stress-energy-momentum nature of the field, we come to the conclusion that Maxwell stress tensor 
M(E, B) should play the basic role, and its zero-divergence in the static case should suggest how 
to determine the appropriate structure and allowed dynamics. 

Recall that any member of the family 



generates the same Maxwell stress tensor. So, the most natural assumption should read like this: 
any member (E, B,Q!i) of this a- family is looking for an energy- momentum exchanging 
partner (E,B,q;2) inside the family, and identifies itself through appropriate (local) 
interaction with the partner found, defining in this way corresponding dynamical 

behavior of the field. 

Simply speaking, a time-dependent EM-field is formally represented by two members of the 
above a-family, and the coupling (E, B,ai) ■<->■ (E,B,q;2) is unique. 

Assuming for simplicity ai = then {£, B, 0) — )■ (E, B) and, making use of the local divergence 
of the Maxwell stress tensor and the time derivative of the local momentum flow of the field, we 
are going to find the value a2 determining the corresponding partner-subsystem. 

Further we shall call these two subsystems just partner-fields. 

Recall the divergence 



As we mentioned, in the static case, i.e. when the vector fields (E, B) do not depend on the time 
"coordinate" ^ = ct, NO propagation of field momentum density P should take place, so, at every 



{£, B) = (E, B, a) = (E cos a — B sin a; E sin a + B cos a), a = a(a;, y, z; t), 




[(rot E) X E Ediv E + (rot B) x B Bdiv B] ^ . 
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point, where (E,B) ^ 0, the stress generated forces must mutually compensate, i.e. the divergence 
VjM'^ should be equal to zero: VjM*-' = 0. In this static case Maxwell vacuum equations 

rotE + — = 0, rotB- — = 0, divE = 0, divB = (J.C.M.) 

OK 

give: rotE = rotB = 0; divE = divB = 0, so, all static solutions to Maxwell equations determine a 
sufficient, but NOT necessary, condition that brings to zero the right hand side of the divergence 
through forcing each of the four vectors there to get zero values. 

In the non-static case, i.e. when ^ 7^ 0; ^ 7^ 0, time change and propagation of field mo- 
mentum density should take place, so, a full mutual compensation of the generated by the Maxwell 
stresses at every spatial point local forces may NOT be possible, which means VjM*-' 7^ in gen- 
eral. These local forces generate time-dependent momentum propagation inside the corresponding 
region. Therefore, if we want to describe this physical process of field energy-momentum density 
time change and spatial propagation we have to introduce explicitly the dependence of the local 
momentum P on (E,B), and to express the flow of the electromagnetic energy-momentum across 
an arbitrary static finite 2-dimensional surface S in two ways: in terms of ViM^^ ^ and in 
terms of the time change of P(E, B), and then to appropriately equilize them. Hence, we have to 
construct the corresponding two differential 2-forms to be integrated on S. 

Note that compare to classical approach where the flows of the very E, B through some 2-surface 
are considered, we consider flows of quantities having direct stress-energy-momentum change sense. 

In terms of = ViM^^ ^ the 2-form that is to be integrated on S is given by reducing 
i^Uo = i^{y^dx^ A dx^ A dx^) = *g{'S) on S, where % denotes the interior product between the 
vector field ^ and the volume 3-form ujo, i.e. the local flow of ^ across Uq, and * denotes the 
euclidean Hodge *. This local force flow i^ujo across S generates changes of the momentum density 
flow across S which should be equal to ^ Jg *g(P{'E,'B)) (recall that t is considered as external 
parameter). We obtain (after restricting *^(P(E, B)) and *g{^) on S): 



di 



[ *5(P(E,B))= / *g{S). 
Js Js 



The explicit expression for P(E,B), paying due respect to J.Poynting, and to J.J.Thomson, 
H.Poincare, M. Abraham, and in view of the huge, a century and a half available experience, has 
to be introduced by the following 

Assumption: The entire field momentum density is given by P := ^E x B . 

According to this Assumption, to the above interpretation of the relation ViM*-' 7^ 0, in view 
of the assumed by us local energy-momentum exchange approach to description of the dynamics of 
the field, in vector field terms and in canonical coordinates on we come to the following vector 
differential equation (the 2-surface S is arbitrary and static) 

^(ExB)=5, C = ct, (*) 

which is equivalent to 

(?B \ / 9E 

rotE — X E + EdivE (rotB - — ) xB + BdivB = 0. 
9U V 9C 
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This last equation we write down in the following equivalent way: 



9B 

rotE + — 1 X E + BdivB 
9^ 



dE , 

rotB - — 1 X B + EdivE 
9^ 



(**) 



The above relation (*) and its explicit forms we consider as mathematical adequate in energy- 
momentum-change terms of the electric-magnetic and magnetic-electric induction phenomena in 

the charge free case. We recall that it is usually assumed these induction phenomena to be described 
in classical electrodynamics by the following well known integral equations 

A. 



d 



I *g{B)\s = — *^(rotE)|5 (the Faraday induction law), 
Js Js 

/ *g(E)\s = / *^(rotB)|5 (the Maxwell displacement current law), 
Js Js 



d^Js JS 

where (■••)|s' means restriction of the corresponding 2-form to the 2-surface S. 

We would like to note that these last Faraday-Maxwell relations have NO direct energy-momentum 
change-propagation (i.e. force flow) nature, so they could not be experimentally verified in a direct 
way. Our feeling is that, in fact, they are stronger than needed. So, on the corresponding solutions 
of these equations we'll be able to write down formally adequate energy-momentum change expres- 
sions, but the correspondence of these expressions with the experiment will crucially depend on 
the nature of these solutions. As is well known, the nature of the free solutions (with no boundary 
conditions) to Maxwell vacuum equations with spatially finite and smooth enough initial conditions 
requires strong time-instability (the Poisson theorem for the D'Alcmbcrt wave equation which each 
component of E and B must necessarily satisfy). And time-stability of time-dependent vacuum 
solutions usually requires spatial infinity (plane waves), which is physically senseless. Making calcu- 
lations with spatially finite parts of these spatially infinite solutions may be practically acceptable, 
but from theoretical viewpoint assuming these equations for basic ones seems not acceptable since 
the relation "time stable physical object - exact free solution" is strongly violated. 

Before to go further we write down the right hand side bracket expression in the following two 
equivalent ways: 

di-B] 



rotB -I- 



rot (-B) -I- 



dE 



X B + (-E)div (-E) 



X (-B) EdivE 



These last two expressions can be considered as obtained from the left hand side of the above 
relation under the substitutions (E, B) (B, — E) and (E, B) — >■ (— B, E) respectively. Hence, the 
subsystem (E, B) chooses as a partner (— B,E), or (B,— E), so, 02 = We conclude that the 
subsystem (E, B, a) will choose as partner-susbsystem (E, B, a -|- ^) or (E, B, a — |). 
We may summarize this nonrelativistic approach as follows: 

A real free field consists of two interacting subsystems (Ei, S2), and each subsystem 
is described by two partner-fields inside tiie a{x,y, z;t)-fa.mily 



{S,B) = (Ecosa — Bsina; Esina -|- Bcosa), 



Si = (^aD'Sai), S2 = (£02,^02)^ giving the same Maxwell stress-energy tensor, so the 
full stress energy tensor is the sum: M(Si,S2) = ^M{T,i) + ^M{T,2)- Each partner- 
field has electric and magnetic constituents, and each partner-field is determined by 
the other through (±f ) - rotation-like transformation. Both partner-fields carry the 
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same energy-momentum : M(Ei) = M(S2)? ^^nd the field propagates in space tlirough 
minimizing the relation +/| ^ 0. The intrinsic dynamics of a free real time-dependent 
field establishes and maintains local energy-momentum exchange partnership between 
the two partner-fields, and since these partner-fields carry always the same stress- 
energy, the allowed exchange is necessarily simultaneous and in equal quantities, so, 
each partner-field conserves its energy-momentum during propagation. 



5 Internal interaction and evolution in energy-momentum terms 



In order to find how much is the locally exchanged energy-momentum we are going to interpret ap- 
propriately the above equation (**). Our object of interest, representing the integrity of a real time 



dependent electromagnetic field, is the couple 



(the other case (E,B);(B,-E) 



IS 



(E,B);(-B,E 

considered analogically). In view of the above considerations our equations should directly describe 
admissible energy-momentum exchange between these recognized two subsystems. Since these two 
subsystems are supposed to keep their recognizabilty during spacetime propagation, we have to 
take in view corresponding admissible, i.e. not leading to nonrecognizability but real, changes. It 
is supposed, of course, that these admissible real changes must be mathematicaly represented by 
tensor objects as it is, for example, how a vector field V changes along itself : if the change of V 
is W, then the corresponding "projection" VyF of V on W measures the admissibility, and if 
VyF = the change is admissible. 

Let's denote such changes by D(E,B) and B,E) for each partner-field. Their admissible 
nature now should formally mean that the values of the corresponding "projections" 

<P : D{E, B) ^ (E, B) and ^ : D{-B, E) ^ (-B, E) 

might be zero, i.e. not essential. Moreover, any available "cross-projections" 

: D{E, B) ^ (-B, E) and ^ : D{-B, E) ^ (E, B) 

must be appropriately equalizable, guaranteing in this way the required recognizability of each 
partner-field. Prom physical point of view the values of all these "projections" should have energy- 
momentum change nature in order to be considered as really available. 

The left hand sides of the vacuum Maxwell equations (J. CM.) we consider now as non-zero 
admissible changes, i.e. their projections on the partner-fields might be zero or correspondingly 
equalizable . 

Following this assumption, the change object L'(E,B) for the first partner-field (E, B) we 
naturally define as 

dB 



D(E,B) := ( rotE-F — ; divB 



The corresponding "projection" *p : D(E, B) (E, B) 
q}[D(E,B);(E,B)] =q3 



(dB \ 
rotE+-^; divBj ;(E,B) 



is suggested by the left hand side of the above energy-momentum exchange expressions (**), so we 
define it by : 

(dB \ 1 / dB \ 

rotE-f — ; divBj ;(E,B) :=frotE + -^j xE + BdivB. 
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For the second partner-field (— B,E), following the same procedure we obtain: 



^[D(-B,E);(-B,E)]=^ 



rot(-B) + — ;divEj ;(-B,E) 



(9E \ / 9E \ 

rot (-B) + — j X (-B) + Ediv E = (rot B - — j x B + Ediv E. 

Hence, relation (**) takes the form 

^ [D{E, B); (E, B)] = [D{-B, E); (-B, E)] . 

We turn now to the mutual energy-momentum exchange between the two partner-fields (E, B) 
(— B,E), keeping in mind that each of the two partner fields must keep its energy-momentum . 
The formal expressions are easy to obtain. In fact, in the case £)(— B,E) — >■ (E, B) we have to 
"project" the change object for the second partner-field given by 

L>(-B,E) := (^rot(-B) + ^; divE 

on the first partner-field (E, B). We obtain: 

OE \ ( 9E \ 

rot(-B) -h — j X E-l-BdivE = - f rotB - — j x E-|-BdivE . 

In the reverse case we have to "project" the change-object for the first partner-field (E, B) given 

by 

L>(E,B) := ^rotE + divB^ 

on the second partner-field (— B,E). We obtain 

5B \ / c?B \ 

rot E j X (-B) Ediv B = - (rot E + — j x B Ediv B. 

We recall now the special kind of dynamical equilibrium between the two partner-fields, namely, 
the two partner-fields necessarily carry always the same energy and momentum: 

M^^(E, B, ai) = M^^(E, B, aa), where |ai - asl = |, 

so, the mutual exchange is simultanious and in equal quantities. 

Hence, in view of the required recognizability of each of the two subsystems formally represented 
by the zero values of the two "sclf-projcctions" , and the local dynamical equilibrium between 
these two subsystems formally represented by the simultanious and the equal up to a sign mutual 
"projections" during any time period, the final equations read: 



rot E — X E -F- Bdiv B = 0, 
(9E \ 

rotB-— xB-FEdivE = 0, 

dB\ ( aE\ 

rotE — X B - EdivB + rotB - ^ x E - BdivE = 0. 
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The third equation fixes, namely, that the exchange of energy-momentum density between the two 
partner-fields is simultaneous and in equal quantities, i.e. a permanent dynamical equilib- 
rium between the two partner-fields holds. 

Note that, this double-field viewpoint and the corresponding mutual energy-momentum ex- 
change described by the last equation are essentially new moments. 

The above equations can be given the following form in terms of differential forms. If g is the 
euclidean metric let's introduce the following notations: 



Then we obtain 



g{E) = ri, g{B) = p, g-\*n) = g-H*P) = *P- 

5(rotE X E) = i{'E)dr], ^(rotB x B) = i{B)dp, 

5(rotE X B) = i{B)dr], ^(rotB x E) = i(E)d/3, 

Ediv(B) = i{*r))d * 13, Bdiv(E) = i(*"/3)d * ??. 

Under these notations and relations the above three framed equations are respectively equivalent 
to: 

■ * A 77 



z(E)dr/ + z(*/3)d * /3 



drj 



z(B)d/3 -I- z(*"j7)d * ?7 = * A /3 
i{B)dr] - i{*r])d * P -\- i{'E)dp - i{*P)d * V = * ( ^ - ^ ^ P ] ■ 



Denoting now the Maxwell stress tensors of E and B correspondingly by M(E) and M(B), and 
introducing a new stress tensor T = M(E) -|- M(B) — M(E -)- B) we obtain 



div 



M(E) + M(B) 



^ dE ^ dB „ 
0, divT = — xE-— xB. 



6 Some properties of the equations and their solutions 

Consider the second equation and replace (E, B) acording to 

(E, B) (aE - 6B, 5E + aE), 
where (a, b) are two constants. After the corresponding computation we obtain 



-\-ab 



rot B - ^] X B + E div E 



dB 



+ b' 



rot E + ^1 X E + B div B 



rotE-l- — X B - EdivB + rotB - — x E - BdivE 



OE 



+ 



0. 



Since the constants (a, b) are arbitrary the other two equations follow. The same property holds 
with respect to any of the three equations. 

Corollary. The system of the three equations is invariant with respect to the transformation 



(E, B) (aE - 6B, 6E aE). 
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Writing down this transformation in the form 



a 


b 


-b 


a 



(E, B) ^ (E', B') = (E, B).a{a, b) = (E, B 

= (aE — 6B, 6E + aB), a = const, b = const, 

we get a "right action" of the matrix a on the solutions. The new solution (E', B') has energy and 
momentum densities equal to the old ones multiplied by (a^ + b^). Hence, the space of all solutions 
factors over the action of the group of matrices of the kind 



a{a, b) 



a 


b 


-b 


a 



in the sense, that the corresponding classes are determined by the value of (a^ + 6^). 

All such matrices with nonzero determinant form a group with respect to the usual matrix 
product. A special property of this group is that it represents the symmetries of the canonical 
complex structure in M^. 

Clearly, all solutions to Maxwell pure field equations are solutions to our nonlinear equations, 
we shall call these solutions linear, and will not further be interested of them, we shall concentrate 
our attention on those solutions of our equations which satisfy the conditions 



^ SB , 
rotE + — /O, 



aE 

rotB- — /O, divE/0, divB/0. 



These solutions we call further nonlinear. 

We consider now some properties of the nonlinear solutions. 

1. Among the nonlinear solutions there are no constant ones. 

2. E.B = 0; This is obvious, no proof is needed. 

3. The following relations are also obvious: 



„ 9B 
rotE + - 



.B = 0; 



„ 5E 
rotB- — 



.E = 0. 



4. It is elementary to see from the last two relations that the classical Poynting energy- 
momentum balance equation follows. 

5. E2 = B2. 

In order to prove this let's take the scalar product of the first equation from the left by B. We 
obtain 



+ B^divB = 0. (1) 

Now, multiplying the second equation from the left by E and having in view E.B = 0, we obtain 

E.i ^rotE + X B i - E^divB = 0. 
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This last relation is equivalent to 



-B.j (^votB + X e| - E^divB = 0. (2) 

Now, summing up (1) and (2), in view of divB 7^ 0, we come to tiic desired relation. 

Properties 2. and 5. say that all nonlinear solutions are null fields, i.e. the two well known 
relativistic invariants /i = — and I2 = 2E.B of the field are zero, and this property leads 
to optimisation of the inequality if + I2 ^ (recall the cigen properties of Maxwell stress tensor), 
which in turn guarantees a{x, y, z; f)-invariance of Ii = I2 = 0. 

6. The helicity property: 

/ dB\ f 5B\ 

B. I rot B - — \ - E. I rot E + — j = B.rotB - E.rotE = 0. 

To prove this property we first multiply (vector product) the third equation from the right by E, 
recall property 2., then multiply (scalar product) from the left by E, recall again E.B = 0, then 
multiply from the right (scalar product) by B and recall property 5. 

Property 6. suggests the following consideration. If V is an arbitrary vector field on then 
the quantity V.rotV is known as local helicity and its integral over the whole (compact) region 
occupied by V is known as integral helicity, or just as helicity of V. Hence, property 6. says that 
the electric and magnetic constituents of a nonlinear solution generate the same helieities. If we 
consider (through the euclidean metric g) the 1-form g(E) and denote by d the exterior derivative 
on M^, then 

^(E) A dg(E) = E.rotE Ady A dz, 

so, the zero helicity says that the 1-form ^(E) defines a completely integrable Pfaff system: 5(E) A 
d^(E) = 0. The nonzero helicity says that each of the 1-forms g(E) and ^(B) defines non-integrable 
1-dimensional Pfaff system, so the nonzero helicity defines corresponding curvature. Therefore the 
equality between the E-helicity and the B-helicity suggests to consider the corresponding integral 
helicity 

~g(B)Ad~g{B) = [ ~g{B) Ad~g{B) 

(when it takes finite nonzero values) as a measure of the spin properties of the solution. 

We specially note that the equality of the local helieities defined by E and B holds also, as it is 
easily seen from the above relation, for the solutions of the linear Maxwell vacuum equations, but 
appropriate solutions giving well defined and time independent integral helieities in this case are 
missing. The next property shows that our nonlinear solutions admit such appropriate solutions 
giving finite constant integral helieities. 



7. Example of nonlinear solution(s): 
E = 



z z 
X, y, ^ ± z)cos(— K— — h const), (pix, y, C ^ •2)sin(— /t— — h const), 



B = 



Co 

z z 
ib(^(x, y,^ziz z) sin(— K— — h const), ^(j){x, y, ^ ± z)cos(— k— — h const), 



where (l){x, y,^^ ziz z) is an arbitrary positive function, < £0 < cx) is an arbitrary positive constant 
with physical dimension of length, and k takes values ±1. Hence, we are allowed to choose the 
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function to have compact 3d-support, and since the energy density of this solution is 4>^dxAdyAdz 
and the total energy is conserved according to property 4., these solutions will describe time-stable 
and space propagating with the speed of light finite field objects carrying finite integral energy. 
Modifying now the corresponding helicity 3-forms to 

A d5(E) = ^^(B) A d~g{B), 

then the corresponding 3d integral gives kTE, where k = ±1, T = 27r Cg/c and E = J cp'^dxAdyAdz 
is the integral energy of the solution. 



7 Scale factor 

We consider the vector fields 

-J ^ dB E X B , . „ 
^ = rotE + - + ^-^d.vB, 

^ dE E X B ^ 
7W = rot B - — - — — — divE, 
|E X B| 

defined by a nonlinear solution. 

It is obvious that on the solutions of Maxwell's vacuum equations T and M are equal to zero. 
Note also that under the transformation (E, B) (— B, E) we get T —A4 and A4 ^ T. 

We shall consider now the relation between ^ and E, and between M. and B on the nonlinear 
solutions of our equations assuming that T ^0 and M. j^O. 

Recalling E.B = we obtain 

(E X B) X E = -E X (E X B) = -[E(E.B) - B(E.E)] = B(e2), 
and since |E x B| = |E||B||sin(E,B)| = E^ = B^, we get 

(9B \ 
rotE + — j X E + BdivB = 0, 

according to our first nonlinear equation. 

In the same way, in accordance with our second nonlinear equation, we get x B = 0. In 
other words, on the nonlinear solutions we obtain that J-" is co- linear to E and Ai is co- linear to 
B. Hence, we can write the relations 

f=fi.E, M = f2.B, 

where /i and /2 are two functions, and of course, the interesting cases are /i ^ 0, oo; /2 ^ 0, oo. 
Note that the physical dimension of /i and /2 is the reciprocal to the dimension of coordinates, i.e. 
[fi] = [f2] = [length]-\ 

Note also that ^ and M are mutually orthogonal: T.M = 0. 

We shall prove now that /i = /2- In fact, making use of the same formula for the double vector 
product, used above, we easily obtain 

/xB + 7WxE = 

(5B\ / SEN 

rotE + —J X B + f rotB - — J x E - EdivB - BdivE = 0, 
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in accordance with our third nonlinear equation. Therefore, 



7"xB + A^xE = 

= /lE X B + /2B X E = (/i - /2)E X B = 0. 

The assertion follows. 

The relation \F\ = \M\ is now obvious. 

Note that the two relations = and T.M. = and the duality correspondence 

{J-,M) —7- {—M,J-) make us consider T and M. as nonlinear analogs of E and B respectively. 
These considerations suggest to introduce the quantity 

£(E.B)= ' ' 



l/il I/2I \T\ \M\ 

which we call scale factor. Note that the physical dimension of C is length. Hence, every nonlinear 
solution defines its own scale factor and, concequently, the nonlinear solutions factorize with respect 
to C. It seems natural to connect the constant Co in the above given family of solutions with the 
so introduced scale factor. Assuming C = Co = const, this could be done in the following way. 

A careful look at the solutions above shows that at a given moment, e.g. t = 0, the finite spatial 
support of the function cp is built of continuous sheaf of nonintersecting helices along the coordinate 
z. Every such helix has a special length parameter b = A/27r giving the straight-line advance along 
the external straight-line axis (the coordinate z in our case) for a unit angle, and A is the z-distance 
between two equivalent points on the same helix. So, we may put A = 2TrCo = const, hence, the 
z-size of the solution may, naturally, be bounded by 27: Co- 

Consider now a nonlinear solution with integral energy E and scale factor Co = const. Since 
this solution shall propagate in space with the speed of light c, we may introduce corresponding 
time period T = 2itCo/c, and define the quntity [} = E.T, having physical dimension of "action". 
The temptation to separate a class of solutions requiring f) to be equal to the Planck constant h 
is great, isn't it, especially if this T can be associated with some helix-like real periodicity during 
propagation?! 



8 Conclusion 

The basic idea we followed in this paper was to pass from linear to nonlinear equations in charge 
free classical electrodynamics through giving explicit physical sense of the equations as local stress- 
energy- momentum relations. This view on the subject naturally oriented our attention to the 
quantities carrying the necessary information - the corresponding stress-energy tensors. The ex- 
isting knowledge about the structure and internal dynamics of free electromagnetic fields made us 
assume the notion for two partner- fi,elds internal structure. These two partner-fields carry equal 
local energy densities, and realize local energy exchange without available interaction energy. More- 
over, they strictly respect each other: the exchange is simultanious and in equal quantities, so, both 
partner-fields keep their identity and recognizability. The corresponding internal dynamical struc- 
ture appropriately unifies translation and rotation through unique space-time propagation with the 
fundamental velocity. All Maxwell solutions are duly respected. The new solutions admit finite 
spatial support, and these spatially finite solutions are of photon-like nature: they are time-stable, 
they demonstate intrincically compatible translational-rotatinal dynamical structure, they carry 
finite energy and intrincically determined integral characteristic f) of action nature through appro- 
priate scale factor Co = const. Their integral energy E satisfies relation of the form identical to 
the Planck formula E.T = f). 
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